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Abstract 



^ ■ One of the most common assumptions in the study of neutron star models and their 

I> . oscillations is that the pressure is isotopic, however there are arguments that this may not 

be correct. Thus in the present paper we make a first step towards studying the nonradial 

0\ [ oscillations of neutron stars with an anisotropic pressure. We adopt the so-called Cowling 

. I approximation where the spacetime metric is kept fixed and the oscillation spectrum for the 

^ ■ first few fluid modes is obtained. The effect of the anisotropy on the frequencies is apparent, 

^ ■ although with the present results it might be hard to distinguish it from the changes in the 



frequencies caused by different equations of state. 
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1 Introduction 

The discovery of gravitational waves is one of the most important goals of the astrophysics 
nowadays. A lot of effort is being devoted to this problem worldwide. Ground based 
experiments [1]- [1] as well as space missions [5] are planned and some of them are expected 
to be able to give results in the near future. Parallel to these efforts the first steps towards 
a third generation gravitational wave telescope (the so-called Einstein telescope) which is 
supposed to have much higher sensitivity are being undertaken [6J. The reason why the 
gravitational waves are so difficult to detect is that they are extremely weak which requires 
detectors with very high sensitivity and also accurate waveforms of the signal emitted from 
the astrophysical objects. 

One of the promising sources of gravitational waves is the oscillations of neutron stars [7] . 
A lot of effort has been spent in studying the gravitational wave emission of these objects 
but there is still many important unanswered questions. The gravitational wave emission 
by neutron stars is ultimately connected to their interior structure. In order to predict 
accurately enough the characteristics of the gravitation waves we need adequate relativistic 
models of the neutron star interior. However, at present little is known about the properties 
and the behaviour of matter at very high densities and pressures. So, in modeling the 
neutron star interior we are forced to make certain assumptions about the properties of the 
neutron star matter. Some of these assumptions seem natural from a physical point of view, 
however, there are always uncertainties and suspicions that the assumptions may be not 
fully correct. As the history of science shows, there are surprises sometimes - Nature does 
not always share our notions of what is "natural" . That is why the alternatives should also 
be investigated. 

One of the widely accepted assumptions in studying the equihbrium configurations of 
neutron stars and their oscillations is that the pressure of the neutron star matter is isotropic. 
There are however arguments that the pressure could be anisotropiciO (8]. Some theoretical 
investigations [9|[T0] show that the nuclear matter may be anisotropic at very high densities 
where the nuclear interactions must be treated relativistically. Anisotropy in the fluid 
pressure can be caused by many other factors. Anisotropy can be yielded by the existence 
of a solid core or by the presence of superfluid [H]- [I3], by pion condensation [H], by 
different kind of phase transitions [TB], by the presence of strong magnetic field [T7] or by 
other factors |8] . From a formal point of view the mixture of two fluids is mathematically 
equivalent to an anisotropic fluid [8j, |18| . 

During the last decades, starting with the pioneering work of [19], there have been 
many papers studying anisotropic spherically symmetric static configurations within general 
relativity ^U\- |13]- These studies show that the anisotropy may have non-negligible effects 
on the neutron star structure and properties. For example the anisotropy may influence 
notably the maximal equilibrium mass, maximum redshift and maximum compactness of 
the stars [20], [35]. It is worth noting also that even for stable configurations the anisotropy 
can support outwardly increasing energy density in the star core [42]. 



Generally speaking the anisotropic fluid has pressures which can differ among the spacial directions. 



The fact that the anisotropy can seriously affect the interior structure, and the properties 
of the stellar configurations make us think that the anisotropy may also have a serious 
influence on the gravitational wave emission and more precisely on the gravitational wave 
spectrum of the stellar configurations. Therefore, in the context of the current efforts to 
detect the gravitational waves, it is important to study the gravitational wave spectrum 
of the anisotropic neutron stars. Such a study is twofold. On the one hand, it can reveal 
the basic characteristics of the gravitational wave spectrum of the anisotropic stars and 
the differences with the spectrum of the isotropic stars. On the other hand, such a study 
provides us with a tool to study the reverse problem - to put constraints on the amount of 
neutron star anisotropy using the observed gravitational wave spectrum in the future. 

In the present paper we undertake the first step towards the study of the gravitational 
wave spectrum of anisotropic neutron stars. More precisely, we investigate the spectrum of 
the nonradial oscillations of anisotropic neutron stars in the Cowling approximation. The 
paper is organized as follows. In section 2 we numerically construct equilibrium configura- 
tions describing anisotropic stars. Section 3 is devoted to the derivation of the perturbation 
equations of the anisotropic neutron stars in the Cowling approximation and the formula- 
tion of the boundary value problem for the oscillation spectrum. In section 4 we present 
the numerical results for the oscillation frequencies. The paper ends with conclusions. 

2 Equilibrium anisotropic configurations of neutron 
stars 

In the spherically symmetric casq^l which we will consider in the present paper, the fluid 
anisotropy means that the radial pressure p differs from the transverse pressure q. The 
mathematical description of an anisotropic fluid in spherical symmetry is given by the 
following energy-momentum tensor 

Tf,u = pu^Uy + pkf,ky + q {gfiu + Uf,Uu - kf^K) , (1) 

where g^j^y is the spacetime metric, u^ is the fluid 4-velocity, p is the fluid energy density 
and k^ is the unit radial vector {k^k'^ = 1) with u^^k^ = 0. Note that Qf^^ + u^u,, — kf^k,y is 
the projection tensor onto the 2-surfaces orthogonal to both m^ and k'^. At the center of 
symmetry the anisotropic pressure must vanish since k^ is not defined there. 

For spherically symmetric spacetimes the metric can be written in the well-known form 

ds^ = -e^^de + e^^rfr^ + r'^dd'' + r^ sin^ ed(l)\ (2) 



■I- The assumption of spherical symmetry is apphcable to the static case of matter sources with an energy- 
momentum tensor satisfying \Tg — Tf | << Tg . So the assumption of spherical symmetry is applicable, for 
example, to the cases when the anisotropy is yielded by the existence of a solid core, by the presence of 
superfluid, or by the presence of pion condensation. Also, the spherical symmetry can be used when the 
anisotropy is yielded by weak enough magnetic field. However, if the magnetic field is very strong, as in the 
case of magnetars, the spherical symmetry assumption may not be a good approximation [17] . 



The Einstein field equations 



then reduce to 



Rfiiy - -Kai^yR = ^TiT^jy (3) 



H^e-- + i(l-e-*)=8.,. (4) 

^e-- - 1 (1 - e--) = 8.P, (6) 

while the contracted Bianchi identity 

V^T^, = (6) 

gives 

p' = -{p + pW-^-^, (7) 

where a = p — q. Introducing the local mass m{r) = | (l — e~^^) and expressing $' from 
(|5]) we can write the dimensionally reduced equations in the Tolman-Oppenheimer-Volkoff 
form 



m' = Anpr'^, (8) 

, , , Anpr^ + m 2a .^, 

P' = -{P + P) f ^ . • 9 

r [r — 2m) r 

In order to close our system we should specify the equations of state for p and for a. 
For the radial pressure we will consider a barotropic equation of state and more precisely 

where po = 1-67 x 10^^ g/cm^ in units where c = 1 and we have chosen F = 2.34 and 
K = 0.0195 obtained when fitting the tabulated data for EOS II [44J. The results are 
qualitatively the same for other values of K and F. 

As explained in [42j we cannot just take a = a{p) because this equation of state is 
too restrictive. Instead we should consider quasilocal equation of state a = a{p, p) where 
/i denotes a quasilocal variable. In principle the equation of state a = a{p, p) should 
be determined by the microscopic theory. Unfortunately, at present we do not have a 
good enough microscopic theory to allow us to find the explicit form of the dependence 



a = cx^p, fi). A drawback of the available microscopic models is the fact that the models are 
developed in fiat spacetime and then the results are transferred to curved spacetime which 
is not completely satisfactory [12]. Within the framework of these types of microscopic 
models it is impossible to find the influence of the curved geometry on the equation of state. 
That is why our approach in the present paper is phenomenological. Following [42j for the 
quasilocal variable we take the local compactness /i = = 1 — e^^^ and we consider the 

following equation of state 

a = Xpfi, (11) 

where A is a parameter. Since the local compactness is zero at the center, this guarantees 
that cr(0) = 0. In order to roughly estimate the range of the parameter A we use the results 
of [T5j where the anisotropy is caused by a pion condensation. In [TJ] it is found that 
< a/p < 1 and therefore we could expect that the maximum value of A is of order of 1. 
In the present paper we adopt the range A G [—2, 2]. 

In order to obtain the background solutions which will be perturbed we solve the reduced 
field equations dl]),© and (jTj) with the appropriate boundary conditions 

A(0) = 0, p(0)=po, $(oo)=0. (12) 

The normalized density p, the radial p, and the anisotropic a pressure as functions of 
the radial coordinate r are shown in Figs. [T] and [2] for several neutron star solutions. The 
central energy density po is the same for all of the solutions presented in the figures and the 
results for several values of the parameter A, which controls the anisotropic pressure, are 
shown. It is interesting to note that for large values of A and for large masses, neutron star 
solutions exist for which the density p is not a monotonic function of the radial coordinate 
but has a maximum and these solutions are dynamically stable |42j . 

In Fig. [3] the mass M of the anisotropic neutron stars is shown as a function of the 
central density po and of the radius R for several values of the parameter A. As we can 
see the properties of the star vary significantly when the anisotropic pressure is varied, i.e. 
when we vary the parameter A. The dynamical stability analysis shows that the solutions 
are stable up to the maximum mass of the sequences 



3 Perturbation equations in the Cowling approxima- 
tion 

In this section we derive the equations describing the nonradial perturbations of the anisotropic 
stars in the so-called Cowling approximation |3S], [IH]. In the Cowling approximation the 
spacetime metric is kept fixed. Despite this simplification the Cowling formalism turns out 
to be accurate enough and reproduces the oscillation spectrum with good accuracy. In fact, 
the comparison of the oscillation frequencies obtained by a fully general relativistic numer- 
ical approach and by the Cowling approximation shows that the discrepancy is less than 
20% for the typical stellar models [17] . 
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Figure 1: The normalized energy density p as a function of the radial coordinate r. 
The results for several neutron-star solutions with the same central energy density 
Po = 7.455 X 10^^ g/cm^ and different values of the parameter A are shown (this central 
energy density gives neutron star with mass M = 1.4Mq in the case with zero anisotropic 
pressure, i.e. when A = 0). 
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Figure 2: The radial p and the anisotropic pressure a as functions of the radial coordinate 
r, normalized to the value of the radial pressure at the center of the star pq. The results 
are for the same solutions as in Fig. [H 




Figure 3: The mass of the neutron stars as a function of the central density (left panel) and 
of the radius (right panel) for several values of the parameter A. 

The equations describing the perturbations in the Cowling formalism are obtained by 
varying the equations for the conservation of the energy-momentum tensor ([6]). Taking into 
account that the metric is kept fixed, we find V u^T^ = where 



5r; = {Sp + Sq) uX + ip + q) iu^.Su'' + SuX) + SgS;: + Sak^k^ + aSk'^k^ + ak'5k^. (13) 
Projecting equation VpST'j^ = along the background 4-velocity m^ we have 

u'VJp + V. { [(p + q)6; + ak'^k^] 6u^} + {p + q) aju" + V,u^5 {ak'k^') = 0. (14) 

Projecting orthogonally to the background 4-velocity by using the operator V^ = ^u + 
ti^M^, we obtain 



{5p + 5q) a^ + {p + q) u" (V,<5m^ - V^<5m,) + V^^g + u^v"^ M + V^^ J (ak'^k,) = 0, (15) 

where a^ = u^V^u^ is the background 4-acceleration. 

At this stage we can express the perturbations of the 4-velocity via the Lagrangian 
displacement vector ^*, namely, 



'dt 



5u' 



where i = 1,2,3 = r,6, (j). 

Now let us consider eq. ( fTSl) ioi p = 9 and p 
we find 



i). Since ae 



(16) 

and uf" = (u\ 0, 0, 0) 



6 



{p + q){uydtie + de5q = Q, 



(17) 



[p + q)[uydtU + d^6q = Q. 



(18) 



Taking into account that p, q and m* depend on r only, the integrabihty condition for the 
above equations gives 



deU = d^ie- 



(19) 



From this condition and the fact that the background is spherically symmetric we find that 
^0 and ^^ are of the form 



Im 

e</, = -$^v^™(r,t)9^y, 



my^ 1 'rji 



(20) 
(21) 



Irn 



where Yim{0, <p) are the spherical harmonics. From now on, in order to simplify the notations, 
we will just write ^ = —VYim when we have expansion in spherical harmonics. 

We proceed further with finding the expressions for the density and pressure perturba- 
tions. From eq.f ll4p after some algebra we find 



1 



-A{V^g[{p + q)C + cr{k^i'W]] 



6p = 

[ip + q)e + crik,e)k'] d, ln(n*) - (p + p)a£ 
It is convenient to express EJ' in the form 



r 



-^y, 



Im 



(22) 



(23) 



and substituting in the above equations, after some algebra we finc^l 



5p = -{p + p) 



-,-'-^^'£±}lv 









Y,„ - ^c-*^i',„ + ^e-*»'y,„ + a 



dr r^ 
where in the last step we have taken into account that a,. = $' and eq. 



^-^^vnu^^) 



''The derivative with respect to the radial coordinate r wiU be denoted by prime or by the standard 
symbol interchangeably. 
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In order to find the perturbation of tfie radial pressure we first use the relation between 
the Eulerian and Lagrangian variations, namely 



5p = Ap — CdrV 
with Ap being the Lagrangian variation. From the equation of state we have 



(25) 



Ap = I Ap = I (5P + i''drP) . (26) 

In this way we obtain the following formula for the perturbation of the radial pressure 






a 



- 2— e-'W - a 



l{l + l) 



V ^ ¥,„,. - 



dp ( 
Sp= -— Up + P) 

For the perturbation of the anisotropic pressure a = a{p,fi) we have 

da = -7:-dp, 
op 

where we have taken into account that 6fi = 0. 

The dynamical equations for W and V follow from eq. flTSl) . namely 



dp _^W 
dr r^ 



^«m-(27) 



(28) 



(P + P) — ^d'^W + dr6p + {6p + 6p)ar + -6a = 0, 
[p + p- a)e-^^d^V -6p + 6a = 0, 



(29) 
(30) 



where 6p are the coefficients in the expansion in the spherical harmonics Yim, i.e 6p = 6pYim- 
From now on we will assume for the perturbation functions a harmonic dependence on 
time, i.e. W{r^t) = W{r)e^^^ and V{r,t) = V{r)e^'^^. Then the above equations become 



,A-2<J| 



-(P + P) — ^00 W + dr6p + {6p + 6p)ar + -6a = 0, 



(p + p- a)e-"^uj'V -6p + 6a = 0. 



(31) 
(32) 



The system (!3T]) - (!32!) can be considerably simplified by combining the equations in an 
appropriate manner. Differentiating equation ( 132|) and adding it to equation ( |3T]) . and also 
using eq.©, we find 



V = 2V<^' 



+ 



a' 



p + p~ a p + p- a 



l-^) P + P ^W 
dp J p + p — a r'^ 

r op 



dp J 



d'^a I &^a 



dp- 



;P + 



dpdp 



P 



(33) 



V. 



This equation together with equation fl3ip solved for W, form a system which is equivalent 
to ( PT]) -( 15^ but much more tractable: 



/ _ dp 
dp 


u' 


a 

+ 

P + P 



p + p- a 
P + P 



d. 



a 



op 



l{l + l)e^V (34) 



V = 2V^' 



-fl + ^)W^ + /(/ + l)eV 
r \ dp 



l-^] P + P ^-W 



dp J p + p — a r"^ 



(35) 



+ 



a 



dp 
-P + p — o- ' \dp 

1 



(T 



P + P - a 



$' 



2 da ( ^ do\ fd'^a , d^o , 



r dp 



dp J \ dp"^ dpdp 



V. 



The boundary condition at the star surface is that the Lagrangian perturbation of the 
radial pressure vanishes 



2p + p-(y [. da 
p + p \ dp 



,-^^V +($' + - 



2 o 



r p + p 



e-^ = 0. 



(36) 



The boundary conditions at the star center can be obtained by examining the behaviour 
in the vicinity of r = 0. For this purpose it is convenient to introduce the new functions W 
and V defined by 

W = Wr^+\ V = VrK (37) 

Then one can show that at r = the following boundary condition is satisfied 



W = -IV. 



(38) 



4 Oscillation spectrum of the anisotropic neutron stars 

The oscillation spectrum of the anisotropic neutron stars in Cowling approximation can 
be obtained by solving the differential equations fl34|l -f l35|) together with the boundary 
conditions (|36|) and f p8|) . We have calculated the frequencies of the /-modes and the higher 
fluid modes pi and p2- All of the presented dependences are shown up to the maximum 
mass for the corresponding parameters where the solutions become unstable 
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Figure 4: The /-mode frequency as a function of the square root of the average density 
\/MjB?. The results for several values of the parameter A are shown. 



An empirical dependence between the /-mode frequencies and the average density was 
found in [181^9] for the case of isotropic neutron stars and it is interesting to see if it 
changes in our case. The /-mode oscillation frequencies as a function of the square root of 
the average density are presented in Fig. IHfor several values of the parameter A. The graph 
shows that the dependence does not change significantly for small values of the average 
density. Only for large values of the average density and for large absolute values of A, the 
deviation from the isotropic neutron stars (i.e. when A = 0) is more significant. But still the 
uncertainties in obtaining the coefficients in the empirical dependence in |18], which come 
from varying the equation of state, are comparable with the deviation due to the anisotropic 
pressure. 

The frequency / and the normalized frequency uj^ B? jM as a function of the mass 
are shown in Figs. [5] and [6] for the / and the pi modes. Depending on the sign of A the 
frequencies can be larger or smaller than in the case of isotropic neutron stars. As we can 
see the frequencies can change considerably when we increase the absolute value of A. Also 
for fixed value of A the differences with the isotropic neutron stars are bigger for larger 
masses because in this case the compactness /i, which enters the EOS for the anisotropic 
pressure (TTTIl . is larger. Therefore for large absolute values of A and for large masses the 
oscillation frequencies can differ significantly from the isotropic neutron star case. 

The normalized frequencies of the /, pi and "pi modes as a function of A are shown 
in Fig. [7] where the mass M = 1.4Mq is the same for all the solutions. As we can see 
the changes in the frequencies as we vary A, are similar for all the modes. This behavior 
is qualitatively different from some of the alternative models of neutron stars where the 
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Figure 5: The frequency / as a function of the mass M (left panel) and the normalized 
frequency a; as a function of M (right panel) for the / mode. The results for several values 
of the parameter A are shown. 
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Figure 6: The results for the pi-mode of the same solutions as shown on figure [51 
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Figure 7: The normalized frequency a; as a function of A for fixed value of the mass M 
1.4Mq. The results for the /, pi and p2 modes are shown. 



oscillations frequencies vary more significantly as a function of the corresponding parameter, 
for the higher fluid modes [50]- [52]. 

It is interesting to compare the effects on the oscillation spectrum caused by varying the 
anisotropic pressure and by changing the equation of state of the radial pressure. As we 
said before the dependence between the /-mode frequencies and the average density does 
not change much when we vary the equation of state. The same is true also when we vary 
the anisotropic pressure. But the normalized frequency a; as a function of the mass changes 
significantly when we vary A and the equations of state. This can be seen on Fig. [8] where 
the results for the /-mode oscillation frequencies are presented for two equations of state 
of the radial pressure and for several values of A. The EOS II is the standard polytropic 
equation of state which we used up to now with F = 2.34 and K = 0.0195. The EOS A is 
again a polytropic equation of state where the coefficients F = 2.46 and K = 0.00936 are 
obtained when fitting the tabulated data for EOS A [53] and EOS II is stiffer than EOS A. 
As we can see the presence of an anisotropic pressure changes the frequencies in a similar 
way as changing the EOS, and more precisely positive values of A lead to frequencies similar 
to a softer EOS, and negative values of A lead to frequencies similar to a stiffer EOS. Thus 
the oscillation spectrum of neutron stars with anisotropic pressure can mimic to a certain 
extent the oscillation spectrum of neutron stars with softer/stiffer equation of state. But 
still if we consider strong anisotropic pressure the frequencies can change a lot which is 
hard to be achieved by the standard nuclear equations of state. Thus observing more than 
one fiuid mode of a neutron star can help us to prove or at least set limits on the possible 
existence of anisotropic pressure in the neutron stars. 
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Figure 8: The normalized frequency cj as a function of M for the / mode. The results 
for different values of A and for two polytropic equations of state of the radial pressure are 
shown - a soft equation of state EOS A and a stiff equation of state EOS II. 



5 Conclusions 

In the present paper we study how the possible existence of an anisotropic pressure inside a 
neutron star can change the oscillation spectrum. As a first step we examine the oscillations 
in the Cowling approximation where the metric is kept fixed and within this approximation 
the perturbation equations for the anisotropic neutron stars are derived. 

The background solution are obtained numerically by solving the reduced field equations 
where the equation of state for the radial pressure is polytropic and we use a quasilocal 
equation of state for the anisotropic pressure |1^ . The properties of the obtained solutions 
can differ significantly from the isotropic neutron stars. 

The oscillation spectrum of the anisotropic neutron stars is obtained when solving the 
perturbation equations with the appropriate boundary condition and the results for the 
/-mode and the higher fiuid modes are obtained. It turns out that the dependence between 
the /-mode frequencies and the average density which was obtained in [18] does not change 
much in the presence of an anisotropic pressure. The effect of the anisotropy is more 
evident on other dependences, for example, the normalized frequency as a function of the 
mass changes considerably for anisotropic stars. Thus the observation of more than one 
fluid mode can serve as a test for the existence of an anisotropic pressure in the neutron 
stars. 

We have also compared the effect on the oscillation spectrum caused by the anisotropic 
pressure and by changing the equation of state of the radial pressure. It turns out that for 
negative values of the anisotropic pressure, i.e. for A < 0, the changes in the frequencies are 
similar to what we will obtain in the case of isotopic neutron star with a stiffer equation 
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of state, and for A > the results are similar to the case of an isotopic neutron star with 
a softer equation of state. A more detailed analysis, for example if we drop the Cowling 
approximation or if we consider rotating solutions, may show more differences between the 
oscillation spectrum of isotropic and anisotropic neutron stars and we plan to make such a 
study in the future. 

It would be also interesting to check how a change in the equation of state of the 
anisotropic pressure influences the results because the choice of the quasilocal EOS fITT]) is 
by no means the only possible one. Up to now, however, little is know about the EOS of the 
anisotropic pressure and that is why further studies on the possible quasilocal equations of 
state and their effect on the stellar structure and oscillations are needed. 
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